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Abstract 

Magnetic-field - temperature phase diagrams of the axial next-nearest- 
neighbor Ising model are calculated within the framework of a Landau-type 
expansion of the free energy derived from molecular-field theory. Good qual- 
itative agreement is found with recently reported results on body-centered- 
tetragonal UPd2Si2- This work is expected to also be relevant for related 
compounds. 
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There exists a large class of compounds with the body-centered-tetragonal (bet) structure 
of the generic formula AT2X2, where A represents either f/ or a rare-earth element, T is a 
transition metal, and X is either Si or Ge.0 The materials of interest here have strong c- 
axial anisotropy, with the possibility of long-range Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interactions. It may thus be expected that the axial next-nearest-neighbor Ising (ANNNI) 
model and its extensions are relevant in cases where in-plane interactions are ferromagnetic .H 
It appears to be an empirical rule that such models yield a succession of principal phases 
characterized by increasing periodicities as the temperature is increased. A counter example 
is found in the compound UN 128121 where a more complicated model was required to explain 
its observed sequence of magnetically ordered states (Ref. |I], hereafter referred to as I). 
Mean-field theory proved inadequate in this instance. All other relevant AT2X2 compounds, 
however, appear to exhibit magnetic structures consistent with expectations based a mean- 
field approximation (see I). In the present work, the magnetic- field - temperature phase 
diagram of the ANNNI model is considered within the framework of a Landau-type free 
energy developed from molecular-field theory. Good qualitative agreement is found with 
recently reported results on [/P(i25'i2-S^ 

Magnetic interactions are assumed to be well described by a Heisenberg-type Hamiltonian 
of the form 

^ = E Ji^^ - r.)s(r.) • s(r,) - H • ^ s(r,), (1) 

ij i 

where the spin density s(rj), as well as the applied magnetic field H, are assumed to lie 
along the bet c-axis (||z). With ferromagnetic in-plane interactions Jo > 0, the problem 
effectively becomes one- dimensional within mean-field theory. It is the magnetic ordering 
between the ferromagnetic planes separated by c' = |c which is of interest. A Landau-type 
free energy functional may be derived from molecular-field theory, based on the analysis of 
Bak and von Boehmi (also see references in I), with results to sixth-order written as 

F = {1/2V) J dzdz'A{T)s{z)s{z') + (B/AV) J dz[s{z)]^ 

+ {C/6V) J dz[s{z)f + H J dzs{z), (2) 
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where r = z — z' and 



A{z)=aTb{z)-3''J{z)lV, 



(3) 



B = bT, C = cT, with the coefficients a, b, c given in terms of the angular momentum j 
through the Brillouin function (see I). For simphcity, we take here the classical limit j ^ oo 



The ANNNI model with antiferromagnetic next-nearest-neighbor coupling, J2 < 0, is 
frustrated. Numerous higher-order commensurate phases, as well an incommensurate state 
(IC), result from mean- field analyses of a model with near-neighbor interactions also anti- 
ferromagnetic, Ji < 0, and with an additional small ferromagnetic third-neighbor coupling, 
J3 > 0.i Monte-Carlo simulations seem to suggest, however, that only a few princicpal com- 
mensurate phases survive the effects of critical fluctuations.Bi For values of J2/J1 not too 
large, only two commensurate states appear: The period-2 (1) state and the period-3 (12) 
state (see Ref. § for an explanation of the notation). It is precisely these two commensurate 
states which occur in the axial AT2X2 compounds (see I). 

With the limited goal of describing these two commensurate, as well as the IC, states, 
the spin density can be written as 



where m is the uniform component, S is the complex polarization vector, and Q is the wave 
vector. Using this expression in (2) yields many types of Umklapp terms, non-zero only 
for nQ = G, where G is a reciprocal lattice vector. For the purposes of the this work, 
contributions with n = 1 and n = 5 are omitted. It is instructive to present the many terms 
which are relevant, and to write these as 



(with energies divided by j^) so that a = 3,b 



|,c~ 1.697. 



s{z) =m + Se'^' + S*e~'^\ 



(4) 



(5) 



with 



F2 



^Aq[S^ + C.C.]A2Q,G, 



(6) 
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F4 = lB{S^ + 2 I ^2 |2 I |2 

+ 2[2{mSy + S^S^ + c.c.]A2Q,G 

+ A[mS^ + c.c.]A3Q,G + [S' + c.c.]A4q,g}, (7) 



Fg = lC{S^ + 6S^ I |2 +24^1 I mS ^ +12[m^S\S*f + c.c] 
+ 3[4SUmSy + S^S^ + 8\mS\^ \ \^ + c.c.]A2q,g 
+ 2[6S^mS^ + A{mSf + 3{mS*)S^ + c.c.]A3q,g 

+ 3[^|5^ + 4m^S' + c.c.]A4Q,G + [S' + c.c.]A,q,g}, (8) 

where = + 2 | S" p. The coefficients in (6) are given by Aq = aT — Jq, where Jq is the 
Fourier transform of the exchange integral, 

Jq = 4 Jo + 2[JiCOsq + J2COs{2q) + J3COs(3g)], (9) 

with q = cQ. The commensurate states (1) and (12) are represented here by g = vr and 
q = Ivr, respectively. The IC phase is determined by the value of q which maximizes Jq. 
This is usually the ffist ordered state which occurs as the temperature is lowered from the 
paramagnetic (P) phase, with the Neel temperature T/vi = Jq/a. Due to the Umklapp 
terms, the free energies of each of these ordered states must be evaluated separately. None 
of the Umklapp terms contribute to F of the JC-phase. All non-Umklapp terms, as well 
as those with n = 2, 4, and 6, must be included in the free energy for the (l)-phase, 
whereas only non-Umklapp terms and those with n = 3 contribute to the (12)-phase. These 
separate free energies must be minimized as a function of the relevant variables in each case: 
Ficijn, I S |,g), F^i)(m, S"), and F(^i2){m, \ S where S* =| 5 | e*"^. The stable phase as 
a function of temperature, magnetic field and exchange parameters, is then determined by 
comparing the numerical results for these free energies. For the (12)-state, the free energy 
is minimized by a phase angle = {2m + 1)tt/3, where m is an integer. 

It is convenient to consider phases which occur as a function of J2/J1, with J3/J1 set to 
a small value.iHl In the absence of a magnetic field, a number of sequences of phases results 

4 



from the present model depending on J2. With increasing temperature, the following can 
occur (see I) (1) ^ P, (1) ^ IC ^ P, (1) ^ (12) ^ IC ^ P, (12) ^ IC ^ P. The 
resulting J2/ Ji — T phase diagram corresponds qualitatively to the Monte-Carlo simulation 
results (see Fig. 1 of Ref. ||). (The sequence (12) (1) IC P observed in UNi2Si2 
cannot be explained by mean-field theory.) All transitions except for IC — P are first order. 
It is important to note that the n = 3 Umklapp terms in (7) and (8) are linear in m. As 
a consequence, the spin density of the (12) phase has a non-zero uniform component even 
in the absence of an applied field. It is this coupling between m and S which enhances the 
stablity of this state in the presence of a magnetic field. 

Magnetic phase diagrams were calculated in the manner described above. Exchange 
parameters were set to Jq = 1, Ji = —1, and J3 = 0.03. With these values, phases (1) and 
(12) are degenerate in energy at T = withi'i J2 <^ 0.35 and the ordering wave vector in the 
IC phase is approximately |7r, as observed experimentally in some AT2X2 compounds (such 
as U Pd2Si2). There is no temperature or field dependence of the IC wave vector within the 
present model. This may be included by adding biquadratic exchange (see I). 

Typical results shown in Figs. 1 and 2 demonstrate the increasing stability of the (12)- 
phase with increasing field. (This feature is also observed in UN 128x2^) The phase diagram 
of Fig. 1 corresponds qualitatively to experimental results on U Pd2Si^. Only relatively 
low field strengths were available in this experiment and the merging of the (12) — IC and 
IC — P boundary lines was not observed, although there is clear indication of this tendency. 
For smaller values of J2I J\i the (12)-phase may not occur at available field strengths. At 
even smaller values of this parameter, the JC-phase also disappears. For larger J2/J1, only 
the (12)-state, and perhaps a small region of the JC-state, will occur in the phase diagram. 
At sufficiently large values of this parameter, a period-4 state is stabilized.i'i 

In conclusion, a simple mean-field model based on a Landau-type free energy derived 
from molecular-field theory has been shown to capture the essential features of the ANNNI 
model in a magnetic field. It describes many of the sequences of magnetically ordered states 
which occur in axial AT2X2 compounds. In particular, the magnetic phase diagram of 
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UP 6,2812 is well described by this model. It is hoped that this work will serve to stimulate 
further experimental investigation of magnetic field effects in this class of materials. 
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FIGURES 

FIG. 1. Magnetic phase diagram with exchange parameters Jq = 1, Ji = — 1, J2 = —0.30, and 
= 0.03. Sohd and broken curves represent first-order and continuous transitions, respectively. 



FIG. 2. As in Fig. 1, with J2 = -0.35. 
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